* Wealth index calculation. Trevor Croft, 12/15/2017. Togo MIS 2017. cd 'C:\Users\21180\ICF\Togo MIS 2017 - DP\wealth'. include file='wealth_assets.sps'. missing values QH113 (99.9). execute. set tvars=both. * Macro definitions to perform some recoding automatically. * Note: Do not use an asterisk in front of a macro name to try to comment it out - it will run any way because of the way the macros are expanded. * Delete the line instead or change the spelling of the macro in the comment. * Do not use the names of macros in comments - they will also be expanded too, and likely will fail. *{Create binary variables based on condition and give label }. define dichotomize( !positional !enclose('(',')') / !positional !tokens(1) / !positional !tokens(1) ) compute !2=0. if (!1) !2=1. variable labels !2 !3. value labels !2 0 'No' 1 'Yes'. formats !2 (f1.0). !enddefine. *(Separate shared toilets from non-shared toilets). define shared_fac( !positional !tokens(1) / !positional !tokens(1) / !positional !tokens(1) ). compute !2 = 0. do if (QH106=1). + if (!1=1) !2 = 1. + compute !1 = 0. end if. variable labels !2 !3. value labels !2 0 'No' 1 'Yes'. formats !2 (f1.0). !enddefine. *{Reset missing values to 'does not have', change 2 code to 0}. define no2zero( !positional !tokens(1) ) if (missing(!1) | !1<>1) !1=0. value labels !1 0 'No' 1 'Yes'. !enddefine. *{Construct Variables}. * Source of drinking water. dichotomize (QH101=11) QH101_11 "Source of drinking water: Robinet dans le logement". dichotomize (QH101=12) QH101_12 "Source of drinking water: Robinet dans la cour/parcelle". dichotomize (QH101=13) QH101_13 "Source of drinking water: Robinet chez le voisin". dichotomize (QH101=14) QH101_14 "Source of drinking water: Robinet public/borne fontaine". dichotomize (QH101=21) QH101_21 "Source of drinking water: Puits à pompe ou forage". dichotomize (QH101=31) QH101_31 "Source of drinking water: Puits protégés". dichotomize (QH101=32) QH101_32 "Source of drinking water: Puits non protégés". dichotomize (QH101=41) QH101_41 "Source of drinking water: Source protégée". dichotomize (QH101=42) QH101_42 "Source of drinking water: Source non protégée". dichotomize (QH101=51) QH101_51 "Source of drinking water: Eau de pluie". dichotomize (QH101=61 or QH101=71) QH101_61 "Source of drinking water: Camion citerne ou Charette avec petite citerne/tonneau". *ichotomize (QH101=71) QH101_71 "Source of drinking water: Charette avec petite citerne/tonneau". dichotomize (QH101=81) QH101_81 "Source of drinking water: Eau de surface (Rivière/barrage/Lac/Maré/Fleuve/Canal/Canal d'irrigation)". dichotomize (QH101=91) QH101_91 "Source of drinking water: Eau en bouteille". dichotomize (QH101=92) QH101_92 "Source of drinking water: Eau en sachet". *ichotomize (QH101=96) QH101_96 "Source of drinking water: Autre". * Type of toilet facility. dichotomize (QH105=11) QH105_11 "Type of toilet facility: Chasse d'eau - à un système d'égout". dichotomize (QH105=12) QH105_12 "Type of toilet facility: Chasse d'eau - à une fosse septique". dichotomize (QH105=13) QH105_13 "Type of toilet facility: Chasse d'eau - fosse d'aisance". dichotomize (QH105=14 or QH105=15) QH105_14 "Type of toilet facility: Chasse d'eau - à quelque chose d'autre". *ichotomize (QH105=15) QH105_15 "Type of toilet facility: Chasse d'eau - ne sait pas". dichotomize (QH105=21) QH105_21 "Type of toilet facility: Fosses/latrines - ventilées ameliorées (VIP)". dichotomize (QH105=22) QH105_22 "Type of toilet facility: Fosses/latrines - avec dalles". dichotomize (QH105=23) QH105_23 "Type of toilet facility: Fosses/latrines - sans dalles/trou ouvert". dichotomize (QH105=31) QH105_31 "Type of toilet facility: Toilettes à compostage". *ichotomize (QH105=41) QH105_41 "Type of toilet facility: Seau/tinette". dichotomize (QH105=51) QH105_51 "Type of toilet facility: Toilettes /latrines suspendues". dichotomize (QH105=61) QH105_61 "Type of toilet facility: Pas de toilette/nature". dichotomize (QH105=96 or QH105=41) QH105_96 "Type of toilet facility: Autre". shared_fac QH105_11 QH105_11_sh "Type of toilet facility: Chasse d'eau - à un système d'égout - shared". shared_fac QH105_12 QH105_12_sh "Type of toilet facility: Chasse d'eau - à une fosse septique - shared". shared_fac QH105_13 QH105_13_sh "Type of toilet facility: Chasse d'eau - fosse d'aisance - shared". shared_fac QH105_14 QH105_14_sh "Type of toilet facility: Chasse d'eau - à quelque chose d'autre - shared". *hared_fac QH105_15 QH105_15_sh "Type of toilet facility: Chasse d'eau - ne sait pas - shared". shared_fac QH105_21 QH105_21_sh "Type of toilet facility: Fosses/latrines - ventilées ameliorées (VIP) - shared". shared_fac QH105_22 QH105_22_sh "Type of toilet facility: Fosses/latrines - avec dalles - shared". shared_fac QH105_23 QH105_23_sh "Type of toilet facility: Fosses/latrines - sans dalles/trou ouvert - shared". shared_fac QH105_31 QH105_31_sh "Type of toilet facility: Toilettes à compostage - shared". *hared_fac QH105_41 QH105_41_sh "Type of toilet facility: Seau/tinette - shared". shared_fac QH105_51 QH105_51_sh "Type of toilet facility: Toilettes /latrines suspendues - shared". shared_fac QH105_96 QH105_96_sh "Type of toilet facility: Autre - shared". * Type of cooking fuel. dichotomize (QH108=1) QH108_1 "Type of cooking fuel: Électricité". dichotomize (QH108=2) QH108_2 "Type of cooking fuel: Gaz propane liquéfié (gpl)". dichotomize (QH108=3) QH108_3 "Type of cooking fuel: Gaz naturel/Gaz butane". dichotomize (QH108=4) QH108_4 "Type of cooking fuel: Biogaz". *ichotomize (QH108=5) QH108_5 "Type of cooking fuel: Kérosène/Pétrole". dichotomize (QH108=6) QH108_6 "Type of cooking fuel: Charbon, lignite". dichotomize (QH108=7) QH108_7 "Type of cooking fuel: Charbon de bois". dichotomize (QH108=8) QH108_8 "Type of cooking fuel: Bois". dichotomize (QH108=9) QH108_9 "Type of cooking fuel: Paille/branchages/herbes". *ichotomize (QH108=10) QH108_10 "Type of cooking fuel: Résidus agricoles". *ichotomize (QH108=11) QH108_11 "Type of cooking fuel: Bouse". dichotomize (QH108=12) QH108_12 "Type of cooking fuel: Sciure". dichotomize (QH108=95) QH108_95 "Type of cooking fuel: Pas de repas préparé dans le ménage". dichotomize (QH108=96 or QH108=5 or QH108=10) QH108_96 "Type of cooking fuel: Autre". * Cows/bulls. if (missing(QH111A) | QH110 <> 1) QH111A=0. missing values QH111A (99). compute QH111A_1 = (QH111A >= 1 & QH111A <= 4). compute QH111A_2 = (QH111A >= 5 & QH111A <= 9). compute QH111A_3 = (QH111A >= 10 & QH111A <= 98). variable labels QH111A_1 'Cows/bulls: 1-4' /QH111A_2 'Cows/bulls: 5-9' /QH111A_3 'Cows/bulls: 10+'. * Horses/donkeys/mules. if (missing(QH111B) | QH110 <> 1) QH111B=0. missing values QH111B (99). compute QH111B_1 = (QH111B >= 1 & QH111B <= 4). compute QH111B_2 = (QH111B >= 5 & QH111B <= 9). compute QH111B_3 = (QH111B >= 10 & QH111B <= 98). variable labels QH111B_1 'Horses/donkeys/mules: 1-4' /QH111B_2 'Horses/donkeys/mules: 5-9' /QH111B_3 'Horses/donkeys/mules: 10+'. * Pigs. if (missing(QH111C) | QH110 <> 1) QH111C=0. missing values QH111C (99). compute QH111C_1 = (QH111C >= 1 & QH111C <= 4). compute QH111C_2 = (QH111C >= 5 & QH111C <= 9). compute QH111C_3 = (QH111C >= 10 & QH111C <= 98). variable labels QH111C_1 'Pigs: 1-4' /QH111C_2 'Pigs: 5-9' /QH111C_3 'Pigs: 10+'. * Goats. if (missing(QH111D) | QH110 <> 1) QH111D=0. missing values QH111D (99). compute QH111D_1 = (QH111D >= 1 & QH111D <= 4). compute QH111D_2 = (QH111D >= 5 & QH111D <= 9). compute QH111D_3 = (QH111D >= 10 & QH111D <= 98). variable labels QH111D_1 'Goats: 1-4' /QH111D_2 'Goats: 5-9' /QH111D_3 'Goats: 10+'. * Sheep. if (missing(QH111E) | QH110 <> 1) QH111E=0. missing values QH111E (99). compute QH111E_1 = (QH111E >= 1 & QH111E <= 4). compute QH111E_2 = (QH111E >= 5 & QH111E <= 9). compute QH111E_3 = (QH111E >= 10 & QH111E <= 98). variable labels QH111E_1 'Sheep: 1-4' /QH111E_2 'Sheep: 5-9' /QH111E_3 'Sheep: 10+'. * Chickens or other poultry. if (missing(QH111F) | QH110 <> 1) QH111F=0. missing values QH111F (99). compute QH111F_1 = (QH111F >= 1 & QH111F <= 9). compute QH111F_2 = (QH111F >= 10 & QH111F <= 29). compute QH111F_3 = (QH111F >= 30 & QH111F <= 98). variable labels QH111F_1 'Chickens or other poultry: 1-9' /QH111F_2 'Chickens or other poultry: 10-29' /QH111F_3 'Chickens or other poultry: 30+'. * Guinea fowl. if (missing(QH111G) | QH110 <> 1) QH111G=0. missing values QH111G (99). compute QH111G_1 = (QH111G >= 1 & QH111G <= 4). compute QH111G_2 = (QH111G >= 5 & QH111G <= 9). compute QH111G_3 = (QH111G >= 10 & QH111G <= 98). variable labels QH111G_1 'Guinea fowl: 1-4' /QH111G_2 'Guinea fowl: 5-9' /QH111G_3 'Guinea fowl: 10+'. * Duck. if (missing(QH111H) | QH110 <> 1) QH111H=0. missing values QH111H (99). compute QH111H_1 = (QH111H >= 1 & QH111H <= 4). compute QH111H_2 = (QH111H >= 5 & QH111H <= 9). compute QH111H_3 = (QH111H >= 10 & QH111H <= 98). variable labels QH111H_1 'Duck: 1-4' /QH111H_2 'Duck: 5-9' /QH111H_3 'Duck: 10+'. * Turkey. if (missing(QH111I) | QH110 <> 1) QH111I=0. missing values QH111I (99). compute QH111I_1 = (QH111I >= 1 & QH111I <= 4). compute QH111I_2 = (QH111I >= 5 & QH111I <= 9). compute QH111I_3 = (QH111I >= 10 & QH111I <= 98). variable labels QH111I_1 'Turkey: 1-4' /QH111I_2 'Turkey: 5-9' /QH111I_3 'Turkey: 10+'. * Hectares for agricultural land. compute landarea = QH113. if (missing(QH113) | QH113 >= 99.8) landarea = 99.9. if (missing(QH112) | QH112 <> 1) landarea=0. missing values landarea (99.9). * Electricity. no2zero QH114A. * Radio. no2zero QH114B. * Television. no2zero QH114C. * Telephone (non-mobile). no2zero QH114D. * Computer. no2zero QH114E. * Refrigerator. no2zero QH114F. * CD/DVD player. no2zero QH114G. * Internet connection. no2zero QH114H. * Satellite dish. no2zero QH114I. * Washing machine. no2zero QH114J. * Air conditioner. no2zero QH114K. * Stove/ Gas stove. no2zero QH114L. * Fan. no2zero QH114M. * Watch. no2zero QH115A. * Mobile telephone. no2zero QH115B. * Bicycle. no2zero QH115C. * Motorcycle or scooter. no2zero QH115D. * Animal-drawn cart. no2zero QH115E. * Car or Truck. no2zero QH115F. * Boat with a motor. no2zero QH115G. * Bank account. no2zero QH116. * Main floor material. dichotomize (QH131=11) QH131_11 "Main floor material: Terre/sable". *ichotomize (QH131=12) QH131_12 "Main floor material: Bouse". *ichotomize (QH131=21) QH131_21 "Main floor material: Planches en bois". *ichotomize (QH131=22) QH131_22 "Main floor material: Palmes/bambou". dichotomize (QH131=31) QH131_31 "Main floor material: Parquet ou bois ciré". dichotomize (QH131=32) QH131_32 "Main floor material: Bandes de vinyle/ asphalte". dichotomize (QH131=33) QH131_33 "Main floor material: Carrelage/carreaux". dichotomize (QH131=34) QH131_34 "Main floor material: Ciment". dichotomize (QH131=35) QH131_35 "Main floor material: Moquette". *ichotomize (QH131=96) QH131_96 "Main floor material: Autre". * Main roof material. *ichotomize (QH132=11) QH132_11 "Main roof material: Pas de toit". dichotomize (QH132=12) QH132_12 "Main roof material: Chaume/palmes/feuilles". dichotomize (QH132=13) QH132_13 "Main roof material: Mottes de terre". dichotomize (QH132=21) QH132_21 "Main roof material: Natte". dichotomize (QH132=22) QH132_22 "Main roof material: Palme/bambou". *ichotomize (QH132=23) QH132_23 "Main roof material: Planches en bois". *ichotomize (QH132=24) QH132_24 "Main roof material: Carton". dichotomize (QH132=31) QH132_31 "Main roof material: Tôle". dichotomize (QH132=32) QH132_32 "Main roof material: Bois". dichotomize (QH132=33) QH132_33 "Main roof material: Zinc/fibre de ciment". dichotomize (QH132=34) QH132_34 "Main roof material: Tuiles". dichotomize (QH132=35) QH132_35 "Main roof material: Ciment". *ichotomize (QH132=36) QH132_36 "Main roof material: Shingles". *ichotomize (QH132=96) QH132_96 "Main roof material: Autre". * Main wall material. *ichotomize (QH133=11) QH133_11 "Main wall material: Pas de mur". dichotomize (QH133=12) QH133_12 "Main wall material: Bambou/cane/palme/tronc". dichotomize (QH133=13) QH133_13 "Main wall material: Terre". dichotomize (QH133=21) QH133_21 "Main wall material: Bambou avec boue". dichotomize (QH133=22) QH133_22 "Main wall material: Pierre avec boue". dichotomize (QH133=23) QH133_23 "Main wall material: Adobe non recouvert". dichotomize (QH133=24) QH133_24 "Main wall material: Contre-plaqué". *ichotomize (QH133=25) QH133_25 "Main wall material: Carton". *ichotomize (QH133=26) QH133_26 "Main wall material: Bois de récupération". dichotomize (QH133=31) QH133_31 "Main wall material: Ciment". dichotomize (QH133=32) QH133_32 "Main wall material: Pierre avec chaux/ciment". dichotomize (QH133=33) QH133_33 "Main wall material: Briques". dichotomize (QH133=34) QH133_34 "Main wall material: Blocs de ciment". dichotomize (QH133=35) QH133_35 "Main wall material: Adobe recouvert". dichotomize (QH133=36) QH133_36 "Main wall material: Planche en bois/shingles". *ichotomize (QH133=96) QH133_96 "Main wall material: Autre". *{Members per sleeping room}. if (hhusual=0) hhusual=hhslept. if (QH109>0) memsleep=trunc(hhusual/QH109). if (QH109=0) memsleep=hhusual. if (missing(QH109) or QH109>=99 or memsleep>=98) memsleep=99. variable labels memsleep 'Number of members per sleeping room'. value labels memsleep 0 'Less than 1 per room'. formats memsleep (f2.0). missing values memsleep (99). * Compute urban and rural variables coded (1/0) for filters later. compute urban=(QHTYPE = 1). compute rural=(QHTYPE = 2). variable labels urban 'Urban' / rural 'Rural'. value labels urban 1 'Urban' / rural 1 'Rural'. formats urban rural (f1.0). execute. * Check on indicator variable creation. frequencies variables=QHTYPE to domestic. * Toilet facility by shared/not shared. crosstabs /tables=QH105 by QH106. * Land area by units - if there are separate units - need to convert them to one unit. *crosstabs /tables=QH113N by QH113U. * removed QH101_71 QH101_96 QH105_15 QH105_41 QH105_15_sh QH105_41_sh QH108_5 QH108_10 QH108_11 QH131_12 QH131_21 QH131_22 QH131_96 QH132_11 QH132_23 QH132_24 QH132_36 QH132_96 QH133_11 QH133_25 QH133_26 QH133_96. frequencies variables= QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_2 QH111I_3 landarea . * Turn off weights before all factor analysis. weight off. * Name the dataset window for the HH data for use later. dataset name assets. save outfile='assets.sav'. *** Test Factor Analysis. filter off. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_2 QH111I_3 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_2 QH111I_3 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /method=correlation. *** Common Factor analysis ***. ** Removed area-specific variables **. ** Agricultural animal variables excluded ** Land area excluded ** Any others ?. *** Common Factor Analysis. filter off. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL com) /method=correlation. *** Urban Factor Analysis. filter off. filter by urban. execute. * Removed QH101_42 QH105_31 QH105_96 QH105_31_sh QH105_51_sh QH108_96 QH132_13 QH132_22 QH133_21 QH133_96 QH111A_3 QH111B_1 QH111I_2. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_51 QH105_61 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_21 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 LAND memsleep QH111A_1 QH111A_2 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_3 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_51 QH105_61 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_21 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 LAND memsleep QH111A_1 QH111A_2 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_3 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL urb) /method=correlation. *** Rural Factor Analysis. filter off. filter by rural. execute. * Removed QH105_11_sh QH105_14_sh QH114J QH133_24. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_12_sh QH105_13_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_2 QH111I_3 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_12_sh QH105_13_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_2 QH111I_3 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL rur) /method=correlation. * Label the created score variables. rename variables (com1 urb1 rur1=comscore urbscore rurscore). variable labels comscore 'Common wealth score' /urbscore 'Urban wealth score' /rurscore 'Rural wealth score'. * Add a variable used for linking later. filter off. string rowtype_ (A8). compute rowtype_ = 'EST'. * Calculate regressions. ** Area=urban. filter off. filter by urban. execute. * Declare a dataset to be written to in the regression. dataset declare urbancorv. * Run regression of comscore with urbscore. regression /missing listwise /statistics coeff outs R anova /criteria=pin(.05) pout(.10) /noorigin /dependent comscore /method=enter urbscore /outfile=corv(urbancorv). * Activate file of output from regression. dataset activate urbancorv. * Drop all rows of output except the coefficients. select if (rowtype_ = 'EST'). execute. * Delete unnecessary variables before merging. delete variables DEPVAR_ VARNAME_. * Rename variables containing the constant and the coefficient. rename variables CONST_=urbconst urbscore=urbcoeff. * Re-activate the main household data. dataset activate assets. * merge the coefficients. match files /file = * /table = urbancorv /by ROWTYPE_. execute. * Calculate regressions. ** Area=rural. filter off. filter by rural. execute. * Declare a dataset to be written to in the regression. dataset declare ruralcorv. * Run regression of comscore with rurscore. regression /missing listwise /statistics coeff outs R anova /criteria=pin(.05) pout(.10) /noorigin /dependent comscore /method=enter rurscore /outfile=corv(ruralcorv). * Activate file of output from regression. dataset activate ruralcorv. * Drop all rows of output except the coefficients. select if (rowtype_ = 'EST'). execute. * Delete unnecessary variables before merging. delete variables DEPVAR_ VARNAME_. * Rename variables containing the constant and the coefficient. rename variables CONST_=rurconst rurscore=rurcoeff. * Re-activate the main household data. dataset activate assets. * merge the coefficients. match files /file = * /table = ruralcorv /by ROWTYPE_. execute. dataset close urbancorv. dataset close ruralcorv. dataset activate assets. filter off. *** Calculate combined wealth score from Urban and Rural Scores. * Use coefficients from urban and rural regressions above. compute combscor=0. print formats combscor (f11.5). write formats combscor (f11.5). ** Urban. if (urban = 1) combscor=urbconst+urbcoeff*urbscore. ** Rural. if (rural = 1) combscor=rurconst+rurcoeff*rurscore. variable labels combscor 'Combined national wealth score'. execute. compute hhwt = QHWEIGHT/1000000. variable labels hhwt 'HH weights'. formats hhwt (f12.6). weight by hhwt. filter off. execute. frequencies variables=combscor comscore /format=notable /ntiles=5 /statistics=minimum maximum stddev mean /histogram normal /order=analysis. filter off. filter by urban. execute. frequencies variables=combscor urbscore /format=notable /ntiles=5 /statistics=minimum maximum stddev mean /histogram normal /order=analysis. filter off. filter by rural. execute. frequencies variables=combscor rurscore /format=notable /ntiles=5 /statistics=minimum maximum stddev mean /histogram normal /order=analysis. filter off. execute. *Calculate quintiles and scores for data file. compute hhmemwt=QHWEIGHT*hhusual/1000000. variable labels hhmemwt 'HH members weighting for index'. formats hhmemwt (f12.6). weight by hhmemwt. filter off. filter by urban. execute. rank variables=urbscore (A) /rank /ntiles (5) /print=yes /ties=mean. filter off. filter by rural. execute. rank variables=rurscore (A) /rank /ntiles (5) /print=yes /ties=mean. filter off. execute. rank variables=combscor (A) /rank /ntiles (5) /print=yes /ties=mean. variable labels ncombsco 'Combined wealth index' /nurbscor 'Urban wealth index' /nrurscor 'Rural wealth index'. value labels ncombsco nurbscor nrurscor 1 'Lowest' 2 'Second' 3 'Middle' 4 'Fourth' 5 'Highest'. *** Check on quintiles. frequencies variables=ncombsco nurbscor nrurscor. weight by hhwt. ctables /table ( QH101_11+QH101_12+QH101_13+QH101_14+QH101_21+QH101_31+QH101_32+QH101_41+QH101_42+QH101_51+QH101_61+QH101_81+QH101_91+QH101_92+ QH105_11+QH105_12+QH105_13+QH105_14+QH105_21+QH105_22+QH105_23+QH105_31+QH105_51+QH105_61+QH105_96+ QH105_11_sh+QH105_12_sh+QH105_13_sh+QH105_14_sh+QH105_21_sh+QH105_22_sh+QH105_23_sh+QH105_31_sh+QH105_51_sh+QH105_96_sh+ QH108_1+QH108_2+QH108_3+QH108_4+QH108_6+QH108_7+QH108_8+QH108_9+QH108_12+QH108_95+QH108_96+ QH114A+QH114B+QH114C+QH114D+QH114E+QH114F+QH114G+QH114H+QH114I+QH114J+QH114K+QH114L+QH114M+QH115A+QH115B+QH115C+QH115D+QH115E+QH115F+QH115G+QH116+ QH131_11+QH131_31+QH131_32+QH131_33+QH131_34+QH131_35+ QH132_12+QH132_13+QH132_21+QH132_22+QH132_31+QH132_32+QH132_33+QH132_34+QH132_35+ QH133_12+QH133_13+QH133_21+QH133_22+QH133_23+QH133_24+QH133_31+QH133_32+QH133_33+QH133_34+QH133_35+QH133_36+ LAND+memsleep+ QH111A_1+QH111A_2+QH111A_3+QH111B_1+ QH111B_3+ QH111C_1+QH111C_2+QH111C_3+QH111D_1+QH111D_2+QH111D_3+ QH111E_1+QH111E_2+QH111E_3+QH111F_1+QH111F_2+QH111F_3+ QH111G_1+QH111G_2+QH111G_3+QH111H_1+QH111H_2+QH111H_3+ QH111I_1+QH111I_2+QH111I_3+landarea )[S][Mean F8.3] by (ncombsco+nurbscor+nrurscor)[C] /slabels visible=no. * Alternative instead of ctables. *means tables= QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_14 QH105_21 QH105_22 QH105_23 QH105_31 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_6 QH108_7 QH108_8 QH108_9 QH108_12 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH116 QH131_11 QH131_31 QH131_32 QH131_33 QH131_34 QH131_35 QH132_12 QH132_13 QH132_21 QH132_22 QH132_31 QH132_32 QH132_33 QH132_34 QH132_35 QH133_12 QH133_13 QH133_21 QH133_22 QH133_23 QH133_24 QH133_31 QH133_32 QH133_33 QH133_34 QH133_35 QH133_36 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_3 QH111C_1 QH111C_2 QH111C_3 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_1 QH111G_2 QH111G_3 QH111H_1 QH111H_2 QH111H_3 QH111I_1 QH111I_2 QH111I_3 landarea by ncombsco nurbscor nrurscor /cells mean count stddev. graph /histogram(normal)=combscor /title= 'Distribution of Households by Wealth Scores'. frequencies variables=combscor /format=notable /ntiles=5 /statistics=stddev minimum maximum semean mean median mode skewness seskew kurtosis sekurt /order=analysis. weight off. filter off. write formats combscor urbscore rurscore (f11.5). * Save final dataset of assets. save outfile='assets.sav'. *** Write out scores file. write outfile='scores.dat' table /QHCLUST QHNUMBER combscor ncombsco urbscore nurbscor rurscore nrurscor. execute. * Produce table 2.6 weighted by de jure population to test. weight by hhmemwt. compute total = 0. variable labels total 'Total'. value labels total 0 ' '. ctables /vlabels variables=QHTYPE QHREGION total Ncombsco display=both /table QHTYPE [C] + QHREGION [C] + total [C] BY Ncombsco [C][rowpct.count '%' F8.1, totals [rowpct.count '%' F8.1, count '' F8.0]] /categories variables=Ncombsco total=yes /titles title='Table 2.6 Wealth Quintiles' '' 'Percent distribution of the de jure population by wealth quintiles, according to residence and region, Togo MIS 2017 ' corner='Residence/region' /slabels visible=no.